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The formation of amplitude modulated and phase modulated assemblies of neurons is observed 
in the brain functional activity. The study of the formation of such structures requires that the 
analysis has to be organized in hierarchical levels, microscopic, mesoscopic, macroscopic, each with 
its characteristic space-time scales and the various forms of energy, electric, chemical, thermal 
produced and used by the brain. In this paper, we discuss the microscopic dynamics underlying 
the mesoscopic and the macroscopic levels and focus our attention on the thermodynamics of the 
non-equilibrium phase transitions. We obtain the time-dependent Ginzburg-Landau equation for 
the non- stationary regime and consider the formation of topologically non-trivial structures such as 
the vortex solution. The power laws observed in functional activities of the brain is also discussed 
and related to coherent states characterizing the many-body dissipative model of brain. 

PACS numbers: ll.10.-z, 87.85.dm, 11.30.Qc 

I. INTRODUCTION 

A major effort is under way worldwide to develop superior forms of machine intelligence by applying basic concepts 
and equations from physics to the interpretation of neurobiological and neuropsychological data. The property that 
most clearly distinguishes biological intelligence from contemporary machine intelligence is the rich contextualization 
of information by brains in the construction of knowledge and meaning. Computers and robots operate on information, 
but they don't know what it means. Books and reprints contain information and display it seriatim in accord with 
Shannon's theory of information, but knowledge of what the information means is solely in the brains of the writers 
and readers. We propose that what differentiates knowledge from information are the innumerable linkages among 
myriad fragments of information, which create the moment-to-moment framework for effective intentional action into 
the world [lj. We find that meaning can best be described objectively as created and carried by great fields of neural 
activity, which subjectively we experience as thoughts and perceptions. The contents of such fields are constructed 
from the fragments of information that are imported by sensory neurons and stored by changes in the synaptic linkages 
among the cortical neurons. It is the dynamics of the populations in each sensory cortex that subsequently organizes 
the microscopic fragments into meaningful knowledge by creating macroscopic vector fields of activity that organize 
hundreds of millions of neurons and trillions of synapses. Our aim in this report is to summarize our observations 
of cortical fields in humans and other animals engaged in intelligent behaviors, and then to describe the dynamics 
of cortex using concepts adapted from non-equilibrium thermodynamics Q and quantum field theory 0, 0|- I n 
particular we focus on adapting the time-dependent Ginzburg-Landau equation so as to describe the construction and 
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FIG. 1: Left: The burst of gamma oscillation illustrates the amplitude modulation of the shared carrier wave. Right: AM 
patterns are compared with and without the conditioned stimuli (CS) present in the inhaled air. The change between trial sets 
illustrates consolidation: 'off-line' learning requiring participation of the genome. 



transmission of the wave functions in vector fields that we observe in brain activities and experience as knowledge. 

A thermodynamic model is needed, because intelligence is extremely energy-intensive. Brains consume free energy 
in many forms (electric, magnetic, chemical, metabolic) at rates ten-fold greater than any other organ. Salient among 
these forms is electric current flowing in closed loops that is carried by ions (not electrons) in water. The flux of 
ions within neurons is essential for their functions: to communicate by transmitting pulses with axons, and to sum 
synaptic potentials with dendrites Q . The extracellular flows of loop currents are revealed by ohmic potentials that 
we record from cortical surfaces (the electrocorticogram, ECoG) and from the scalp (the electroencephalogram, EEG). 
The cortex is a thin sheet of neurons covering the outer surface of the brain. The average thickness in humans is 
~ 3 mm, and the surface area is ~ 2,000 cm 2 , a ratio of 1 : 10 5 . Whereas the dynamics of microscopic neurons in 
networks and the scalar potential fields of the EEG/ECoG are described in 3-D, the macroscopic fields of activity 
from 10 8 to 10 9 neurons each supporting 10 4 synapses [6] topologically operate in 2-D, owing to the long correlation 
distances compared with the dimensions of neurons. 

Our standpoint is that cortical dynamics cannot be reduced to the level of neural networks. A macroscopic level of 
mass action emerges from microscopic dynamics, into which the senses inject their microscopic information and from 
which macroscopic neural commands are generated. The ECoG and EEG are appropriate to the macroscopic level, 
because they are summed contributions from interactive masses of neurons. The interactions create dynamic neural 
assemblies with properties that are undetected with pulses from microelectrodes and inexplicable in terms of neural 
nets Q. In the simplest description we conceive cortex as a self- regulating, self-stabilized population of neurons. 
It modulates and is modulated by other parts of the brain, but it does so on its own terms. Whereas the cortical 
neural network dynamics involves pulse and wave frequencies, macroscopic population dynamics uses pulse and wave 
densities, which constitute the state variables in space-time. 

In our analysis it is important that we distinguish the level of function by properly defining the state variables. To 
this aim, since in this paper we are mainly interested in the microscopic dynamics out of which macroscopic structures 
emerge, we consider the many-body dynamics of the vibrational quanta of the electric dipoles of water molecules. 
Water is the matrix in which neurons, glial cells and the whole net of dendrites and axons are embedded, including 
in particular the ionic currents that exercise the electric forces on and in nerve membranes that constitute neural 
activity. The dipole dynamics provides the possibility for long-range correlation, which in turn allows global synaptic 
communication among the neurons, each with every other as required for the construction of meaning. It should be 
clearly stated that the water molecules and their dipole oscillations are not the agency of communication among the 
neurons. As said above, the memory required for macroscopic structure formation is based in the trillions of synapses 
between neurons. Likewise, in our prior publications we have made it clear that the extracellular electric fields of the 
EEG and ECoG are epiphenomenal. They are scalar fields of passive dissipation of electric energy as heat, whereas 
the fields of neural activity are vector fields, which actively derive free energy from chemical energy in every square 
micron and dissipate that energy in clouds of action potentials and ultimately as heat. Our many-body analysis shows 
how such a vector field is indeed generated from the basic dynamics of the dipole vibrational quanta 0, H[ ■ 

The main sources of our experimental data are high-resolution images of the EEG and ECoG from 8x8 planar 
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FIG. 2: A. An example is shown from rabbit ECoG of a cone fitted to the surface given by an 8 x 8 array of values of analytic 
phase with phase lead at the apex ("explosion"). B. Analytic phase from human ECoG with phase lag at the apex ("implosion"). 
C. Summary diagram of the variants of phase patterns seen in cinematic displays of the filtered ECoG. Examples can be seen 
at the URL: http://soma.berkeley.edu 



arrays of 64 electrodes The 64 signals reveal brief epochs (~ 0.1 s) of 3 — 5 cycles of narrow band oscillation. 

The epochs recur 3 — 10 times each second, each time with a different carrier frequency randomly centered in the 
beta-gamma range. The amplitude of the carrier wave is modulated in a spatial amplitude modulated (AM) pattern 
(Fig. 1). The set of 64 amplitudes constitutes a 64 x 1 vector that specifies a point in 64-space. Similar patterns 
accruing from repeated presentation of a conditioned stimulus (CS) form a cluster of points. Differing CSs form 
multiple clusters, one for each CS that a subject can discriminate. The AM patterns lack invariance with respect to 
the CS. The same CS in a different context gives a distinctly different cluster. This contextual dependence shows 
that the AM patterns are dependent on the meaning of the CSs; they are not representations of the CSs but of 
the knowledge the subject has about the CSs. The AM patterns evolve slowly with repetition over time and the 
accumulation of new experiences. Each retrieval of a memory changes the memory by adding new context. 

Each AM pattern is accompanied by a spatial pattern of phase modulation (PM) of the carrier frequency, which 
has the form of a cone (Fig. [2J (see also Appendix [A}. The location and sign of the apex vary randomly from each 
epoch to the next. The gradient of the cone in rad/m is fixed in each epoch and varies between bursts inversely with 
the carrier frequency in rad/s. The ratio gives the phase velocity in m/s, which is determined by the conduction 
velocity of pulses on axons running parallel to the cortical surface. We infer that each new AM pattern forms by 
a state transition, which begins at a site of nucleation on the surface and spreads radially as in the formation of a 
raindrop. The diameter of the AM pattern is determined by the distance at which phase dispersion by the phase 
gradient reduces interaction strength below the half-power level. During an epoch the dynamics is stationary in 
carrier frequency, AM and PM, and the dynamics is near linear. We introduce the Hilbert transform, which gives the 
high temporal resolution needed to display the transitions between successive images [H, [l(| • We also introduce the 
concept of criticality, which we need to help explain the capacity of cortex to undergo the dramatic transitions that 
are required to construct and destroy the images enabling perception Cinematic display of the spatial AM/PM 
patterns reveals repetitive pulsatile expansions or contractions of the wave functions, and in some instances rotation 
either clockwise or counterclockwise indicating the stabilization of the AM/PM patterns by formation of vortices, 
hence vector fields 

The point is that the many-body approach allows the formation of spatially delimited domains that constitute the 
basic facilitating environment wherein the formation occurs of the above described AM and PM structures in the 
cortex. These domains are also temporally limited in time, in the sense that they have a finite lifetime; they condense 
and evaporate continually as in a fog, with power-law distributions of duration and diameter. We treat them as 
the result of a space-time non-homogeneous boson condensation process [l4T - fl7j and as neural avalanches by which 
criticality is maintained |l8|. 

Intermittently we see the long-lasting epochs generated by the long-range correlation of neurons locked in coherent 
oscillation mode carrying AM patterns relating to perception. The dipole quantum field correlation thus acts to 
facilitate expansion of the trillions of microscopic synapses modified by learning into the spatially textured macroscopic 
cognition. The expansion is directed by a mesoscopic Hebbian assembly (Appendix A) into the basin of a macroscopic 
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chaotic attractor that characterizes the AM pattern. The Hebbian assembly provides a burst of transition energy 
that brains require in order to restrict the expansion to intended, learned sensory input and to exclude environmental 
noise. From such a perspective, in the present paper we therefore study the thermodynamics of the phase transition 
processes at the level of the many-body dynamics of the dipole wave quanta with defined topological properties. The 
AM/PM patterns then can be considered to be the macroscopic manifestations of a boson condensation and may 
be characterized by its underlying topological properties. The basic symmetry is the (non-Abelian) rotational SU(2) 
dipole symmetry, and the emerging structure formation is entropic, always incurring loss of information particularly 
in abstraction and inductive generalization. The topology of symmetry breaking from the search/receiving mode to 
the constructive/transmitting mode is related with the singularity in the null spike where prior structure is expunged 
and disorder is initiated, as at center (the core) of a vortex [i[ . 

In Section II we present basic elements of the many-body dissipative model with specific attention to the mechanism 
of boson condensation in a stationary regime. Further details of the formalism of the dissipative many-body model 
can be found in Refs. 0, [H], [H|]- The non-stationary regime is discussed in Section III where we obtain the time- 
dependent Ginzburg-Landau equation. The vortex solution is considered in Section IV and the power laws observed 
in functional activities of the brain are discussed in Section V and related to coherent states characterizing the many- 
body dissipative model. Finally, Section VI is devoted to conclusions. Further remarks on the relevant aspects of the 
observed brain activity and some mathematical details are presented in the Appendices. 



II. THE DISSIPATIVE MANY-BODY MODEL OF BRAIN 

In the dissipative many-body model of brain [Til [l5j , one considers the spontaneous breakdown of the rotational 
symmetry of the dipoles of water molecules. In quantum field theory (QFT) spontaneous breakdown of symmetry 
produces the formation of ordered patterns through the mechanism of boson condensation. Such ordered patterns, 
resulting into "in phase" dipole oscillations and often also space ordering are the manifestation of coherent dynamic 
features, out of which collective macroscopic behaviors of coherent spatial domains appear [2(|. The lagrangian 
describing the dynamics of the system is assumed to be invariant under the SU(2) group of rotations for the molecular 
dipole vibrational field 4>(x,t) = ((f>^(x,t) (j)i(x,t)) T (j =t,4-) denotes the dipole quantum number and assumes, for 
simplicity, only two values "up" and "down") [2lll22|, On the contrary, the ground state (the "vacuum") of the system 
cannot be SU(2) invariant since it has to represent the electret state of water. The SU(2) spontaneous breakdown is 
expressed by the relation 

(0\D^(r,t)\0)=V(v,t)^O, (1) 

where D^(r,t) denotes the dipole density operator in the third direction and P(r,t) is the polarization density. The 
dipole density operators flW(r,t), i — 1,2,3, are functionals of the molecular vibrational field <p{x,t). For example, 
they can be expressed as £)W (r, t) — 4>^(x, t)^<Ji(j)(x, t), i = 1,2,3, with a the Pauli matrices. The conclusions we will 
reach, however, are not dependent on the particular expression of Z)W(r,t) i n terms of <j>(x,t). A general theorem in 
quantum field theory (QFT) [l9|, [2(| I23T - I23 ] shows that spontaneous breakdown of symmetry implies the existence of 
gapless fields (the Nambu-Goldstone modes or particles (NG)). They appear as collective excitations describing the 
long-range correlation among the dipoles generating the ordering of the system (electret state). In the dissipative 
many-body model of brain, a number of functional features of brain can be understood as manifestations of the NG 
collective modes [13, EH ■ 

The spontaneous breakdown of SU(2) symmetry represented by Eq. |T]) means that the vacuum is not invariant 
under rotations induced by D^'(r,t) and D^ 2 \r,t). It is, however, invariant under rotations induced by D^ 3 \r,t), 
i.e. under the U(l) rotation group around the third direction: 

^->cxp(igA 3 (r,t)(T3/2)^ . (2) 

Here q is the electric charge at the head of the molecular dipole. We assume Aa(r, t) — > for \t\ — > oo and/or |r| — > oo. 

Of course, invariance of the lagrangian under local (i.e. with space-time-dependent As(r,t)) U(l) rotation group 
signals the presence (requires the introduction) of the electromagnetic interaction among the dipoles, i.e. the electro- 
magnetic vector potential Afj,(r,t), which transforms as 

A M (r,i) ->A M (r,t)-^A 3 (r,t), (3) 

when the phase transformation around the third axis is performed. In the following, we will use the Coulomb gauge 
V • A = as gauge condition. In terms of the dipole operators, the local U(l) transformation is 



D (±) (r,i) -> exp(Ti<zA 3 M))£> (±) , 



(4) 
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where 

D (± \r,t) ^D [1 \r,t)±iD^{v,t). (5) 

Next we observe that actually the ground state of the system does not exhibit invariance under global (i.e. space- 
time independent A3) U(l) symmetry; this would require the possibility of changing the phase of the molecular 
vibrational field 4>(r,t) simultaneously at every space-point by the constant amount A3. Therefore, the global U(l) 
symmetry is also spontaneously broken. The condition which expresses the symmetry breakdown is the non- vanishing 
expectation value of D^ + \r,t) in the ground state |0): 

(0\D(+\r,t)\0)=v(r,t)^O, (6) 

and its complex conjugate relation involving D^~\r, t). v(r, t) is a complex function and is called the order parameter, 
which is thus a (classical) vector field. It describes macroscopic collective properties of the system. Its space-time 
dependence denotes (space-time) non-homogeneities in the ground state. Since the vector potential is also involved 
in the dynamics, general theorems of QFT predict some dynamical effects globally named as the Anderson-Higgs- 
Kibble mechanism [2(| [27J . Here, we will not consider such a mechanism (see Refs. [lfj HH UHl for a general 
discussion). 

Let P(r, t) and P'(r, t) denote the annihilation and the creation operator, respectively, of the dipole-wave quantum 
(the NG mode implied by the spontaneous breakdown of SU(2) symmetry). By resorting to the results of Refs. 
[HHHI], we have that the ground state turns out to be eigenvector of the dipole-wave operator P(r,t): 

which shows that the ground state is a coherent state: dipole-wave quanta (the NG boson modes) are coherently 
condensate in the state |0) [2l|, |22j. We also have 

(0|I>(-J(r,t)DW(r,t)|0) = 2V{r,t) (0\P\r,t)P(r,t)\0) = \v(r,t)\ 2 , (8) 

We may also define the polarization density as 

(0\D^(r,t)D( + \r,t)\0)=p(r,t)6, (9) 

with p(r, t) the charge density and S the (average) dipole length. By writing the NG condensate density as 
(0|Pt(r,i)P(r,i)|0) = n(r,i), we have p(r,t)S= \v{r,t)\ 2 = 2V(r,t)n(r,t). 
We then write the charge density wave function a(r,t) as 



.(r.^V^Me'^, (10) 

with real p(r,t) and 9(r,t). We have |<r(r,t)| 2 = p(r,t) = N(r,t) in the system ground state. Thus, N(r,t) denotes 
the charge density condensate. Sometimes a(r,t) is also called the macroscopic wave function. We observe that the 
order parameter v(r,t) provides a measure of the NG mode density n(r,t) in the ground state. 
In terms of a, the local gauge transformation (j4]) is 

a{r,t) -> e- iqX3( - r ^a{r,t) , (11) 

Note that the transformation ([J} is induced by the phase transformation 

6(r,t)^9(r,t)-q\ 3 (r,t), (12) 

when Eq. (flQl) is used. One may show [2l| that the phase 0(r,t) represents the NG wave dipole field, also called the 
phason field[20Ll2l|. 

As mentioned, the meaning of the macroscopic wave function, or the order parameter v(r,t), resides in the fact that 
in the presence of the spontaneous breakdown of symmetry, the microscopic quantum components of the system behave 
in collective or coherent way, namely they undergo in phase motion, so that the system is then globally described by 
such a common phase of its quantum constituents. This establishes the link between the system macroscopic wave 
function a(r,t) and the microscopic one, ip(r,t) = \J a(r, t) exp(iS/H), for the quantum components. We assume that 
ip(r, t) satisfies the Schrodinger equation. Thus we will set the phase 0(r, t) = S/H. Then for the momentum P = tov 



6 



we have P = V5 = KV9. In the presence of the A field, the canonical momentum is known to be rnv = P — qA (the 
minimal coupling), thus (cf. Eq. (|B4[0 

_™ V = A --V0. (13) 

q q 

We note that the r.h.s. of Eq. (p~3|) can be considered as a gauge transformation with 6*(r, t) being the gauge function: 
A —> A' = A — ~V#. The requirement that the gauge condition V • A = be invariant under gauge transformations, 

i.e. V - A' =0, implies that 6*(r, t) (and Aa(r, t)) has to be a solution of the equation V 2 # = (and V 2 As(r, t) = 0). In 
such a case, Eq. (|13|) gives V • v = 0. We remark that in general one may consider the so-called boson transformation 

0(r,t)->0(r,f)-c f(v,t), (14) 

with c a convenient constant and f(r,t), called boson transformation function, solution of the equation for 9, 
V 2 /(r,i) = 0, so that Eq. (|T4|) is an invariant transformation for the theory. Such a transformation describes 
the space-time-dependent (i.e. non-homogeneous) boson transformation and it can be shown [1, [HI US US US that in 
order to have observable effects f(r,t) has to carry a topological singularity. This theoretical frame thus predicts the 
appearance of topological structure formation, such as vortices, as further discussed in Section llVI (and in Ref. [4]). 
In the case of the U(l) phase symmetry mentioned above, the stationary function f(x) may carry indeed a vortex 
singularity given by 

f(x) = arctan • 0-5) 

Eq. ([T5jl shows that the phase is undefined on the line r = 0, with r 2 = x\ + x\ (which, as observed in [4], reflects in 
the observed phase indeterminacy in the process of transition between two AM pattern frames). For an analysis of 
the vortex properties associated to Eq. (fT5)) sec Refs. [IS US and the discussion in Section ITVl 

In the following we study the variation in space and in time of the macroscopic wave function and of the phason 
field condensate. In QFT these variations denote transformations through physically different vacua, i.e. unitarily 
inequivalent vacua, or, in other words, they denote "phase transitions" . Of course, going through phase transitions, 
the system moves toward the equilibrium, or stationary regime, which is the one minimizing the free energy density 
functional F(a, a*, A) with respect to a*: dF/da* = 0. However, during the phase transition processes such a 
minimization condition is not satisfied, dF/da* =/= 0. Here we are interested in such dynamical transition processes 
since the brain is a far from the equilibrium system, indeed. 

Let us start by considering the Schrodinger equation for a(r,t): 

ih^- = Ha, (16) 
dt x ' 

with 

H=^- {-ihV - qA) 2 + U, (17) 
2m 

where the potential energy U is 

U = qcp + n . (18) 

Here ip = V/q, with V the scalar potential and /iq the chemical potential. (1q is included in the Schrodinger's equation 
for a in order to account for the variations in the number \a\ 2 = N. In the following our discussion we will closely 
follow Ref. |3lj ]. where the time-dependent Ginzburg-Landau equation for the modulus of the order parameter is 
obtained and conservation laws as well as dissipative effects are taken into account. 

The order parameter and the macroscopic wave function a are not affected by quantum fluctuations (in this sense 
they are macroscopic fields). However, as said above, we are interested in their variations occurring in the (non- 
equilibrium) phase transition processes. Therefore, we need to allow variations in the number N, which justify the 
introduction of the chemical potential [iq in Eq. (fl7|) . We note that a chemical potential term is not included in the 
Schrodinger's equation for the microscopic wave function ^(r, t) describing the elementary dipole components, whose 
number is assumed to be constant during the system time evolution. 

Use of Eq. (TlT)|) into Eq. ([To]) gives the the equations for the imaginary and the real parts: 

dN 

°_ + V • (JVv) = 0, (19) 
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mv dO 

fi +H 1 + q<p+ — + h—=0, (20) 

respectively. Here and in the following v 2 — v • v. Eq. (I19|) is the continuity equation (to be compared with Eq. (|B1[) V 
The additional contribution [i\ to the chemical potential in Eq. (|20[) is 

^^^VF" (21) 
and is due to non- homogeneous density of the condensate. From Eq. (|20|) one obtains (see Appendix |5J 

= g(E + vxB)-V(/i + Aii), (22) 

at 

with 

C7.A 77? 

E = -— -Vy? and B = VxA=--(Vxv). (23) 
We observe that, provided that V • v = (see the remarks after Eq. (fT3"]) and Eq. (|D6[) ). Eq. (|T5)) gives 

which expresses the charge density conservation in the stationary regime. Since TV is related to the order parameter, 
this means that the long range correlation, namely the ordering, is preserved in time. In the frame of the dissipative 
many-body model, such a case describes stationary neuronal correlates. However, observations of the functional 
activity of the brain show that a succession of neuronal correlates ( "wave packets" ) , modulated in amplitudes and 
in phase, occurs. In its functional activity the brain appears undergoing a continuous stream of phase transitions. 
We are therefore interested in studying the non-stationary dissipative phase transition processes, where, contrarily 
to Eq. ([M]) . dN/dt ^ 0. Thus we need to study the non-stationary or time-dependent Ginzburg-Landau equation in 
order to consider the possibility of dissipative processes. 



III. TIME-DEPENDENT GINZBURG-LANDAU EQUATION 



Let us start by considering the Ginzburg-Landau (GL) functional F(a, a* , A) representing the free energy density. 
For our task we do not need to specify the explicit form of F(a, a* , A). This will depend on the particular dynamical 
model one adopts for the description of the system under study. In general, it is a non-linear functional of the fields, 
containing a kinetic energy term (— ifiV — qA) 2 /2m and some potential term. The stationary Ginzburg-Landau (GL) 
equation is obtained by extremizing the free energy density F (a, <r*, A) with respect to a*: 



where, in full generality, we may write: 

OF 
da* 



2m 



dF -0 

da* 



(-iHV - qA) 2 + a 



(25) 



(26) 



with a(T) and (3(T) = \a(T)\/N°(T) > acting as a mass term and a positive coupling constant term, respectively. 
N°(T) = \a\ 2 is the equilibrium condensate density. Their values depend on the temperature T and their explicit 
expressions depend on the system under consideration and on the adopted phenomenological model (in solid state 
physics |a(T)| is related to the chemical potential /j,o, expressing, e.g. in superconductivity [3l[, the vacuum perme- 
ability, and f3(T) to the critical (magnetic) fields acting upon the system. Notice that the contribution (a + f3\a\ 2 )a 
in the r.h.s. of Eq. (|26|) may be considered as derived from the potential V(a,a*) = —(/3 /2)(\a\ 2 + a/(3) 2 . The (mean 
value of the) order parameter a minimizing the potential V(a, a*) is zero (disordered or symmetric ground state) for 
a > 0, and non-zero for a < 0, with \a\ 2 = —a//3 ^ (the ordered or asymmetric ground state) 1 



We also note the strong analogy with the laser phase transition [32l |. where the potential V(cr, a*) is again considered. In the phase 
transition between the disordered and the ordered (laser) state (the lasering process), the order parameter \a\ 2 changes in time from 
zero to a non-vanishing value proportional to a, going through the threshold set at a = 0. In the Haken interpretation, a is the pump 
parameter whose tuning may carry the system in the lasering region; a denotes the amplitude of the classical electromagnetic mode. 
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Eq. (12 5 P expresses the condition for the stationary (equilibrium) dynamical regime, where the evolution of the 
macroscopic wave function is controlled by the Schodinger equation (|16l) . Instead, in the case of non-stationary 
regime, the non-vanishing quantity dF/da* expresses the rate at which cr(r, t) approaches its stationary value at the 
minimum of the free energy. It is then customary to consider the generalized time-dependent GL (TDGL) equation: 

^da - i dF . . 

ih—=Ha--—, 27 
at 7 da* 

where now dF/da* is not zero, 7 is a relaxation parameter and the stationary regime is reached in the limit 
(l/j)dF/da* —> 0. In such a limit, the equilibrium wave function is u°(t) — \o°\ exp(— iet/h), with Ha° — ea°. 
The second term in the r.h.s. of Eq. (|27|) thus describes the dissipative contribution coming from incoherent relax- 
ation processes [3l|. Using a = y/N exp(i6), the imaginary and real part of Eq. (|27|) give (see the Appendix [C|l the 
continuity equation: 

™+V-(Nv) = -^N (28) 

Ot T GL 

and the relation for the chemical potentials 



respectively, where 



mv 2 86 , , 



2~/N t gl q N 



is related to the relaxation parameter 7, is of dissipative nature and proportional to non- homogeneities of the con- 
densate. We have used the definition J = qNv. In Eq. (j2"5)) G is given by 



„ ^ 2m f mv 2 \ N . . 

G = ^ W {— +M 1 J+^-l. (3D 

and tgl and £ GL are the GL relaxation time and coherent length, respectively (see Eqs. (|Cip ). By operating with 
the gradient operator on the Eq. (|2"5)) we obtain 

=?(E + vxB)- V(mo+j"i+j"2). (32) 

at 

Eqs. (|2"5)l and (|3"2")l have to be compared with Eqs. (JTHJ) and (|2"2"j) . respectively. 

By using Eq. (|3~Tj) and V • v = 0, the continuity equation (j2"5)l can be rewritten as 

_,(- + vv)iV = -r R iV^0 (33) 

where we have set 

1 2G 

r R = — = — (34) 

TR TQL 

We thus see that the relaxation term Tvldiss = ^rN accounts for the rate of change of the condensate dN/dt. From 
the expression of G, Eq. (|3"Tj) . we see that we may write F R = Ti + T 2 , where 

_ fp T 2m , „ ^ fmv\ 2 ,„_, 
r i = — T? mv =2D GL — , (35) 



TGL 



h 2 ^ \ h 



with the diffusion coefficient Z?gl = £gl/ t gl, and T2 related with non-homogeneities (/ii and N ^ N°). T2 is related 
with dissipative processes whose life-time is usually longer than the GL relaxation time tgl ■ In general, the condition 
tr = tql/2G tgl, requiring, in order to hold, small values of G, ensures fast formation of quasi-local equilibrium 
in the condensate, controlled by tql, compared with longer decay time of the condensate density N. At a critical 
temperature T c , a fast transition to the equilibrium condensate regime occurs; one may then expect G — ¥ 0, so that 
tgl < -tr = t gl /2G, as T ->■ T c . 
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Eq. ((33]) gives the general TDGL equation for the normalized wave function \ — I "!/! " ! = (N/N°) 



1/2. 



dx 
dt 



1 

TGL 



5GL 



v h 



(36) 



where dx/dt — dx/dt + v • V%. In the approximation of Fr w T\ = 2/te, and in the case of fast formation of 
quasi-local equilibrium, we have 



— 4r T 1 



In such a limit, for tql <C te, Eq. (|36[) becomes 



GL 



V h 



< l. 



„ _2 <^X 2\ X 

■DglV x - — = -(1 - x ) 

dt TGL 



(37) 



(38) 



which further reduces to 



CglV 2 x + X ~ X 3 = 



(39) 



when the stationary condition dx/dt = is met, i.e. dx/dt = and v • Vx = 0. Eq. (|39)) reproduces indeed the 
stationary GL equation (|C6I) when the disequality (|3"T1) holds. 



IV. THE VORTEX SOLUTION 



Our discussion in the previous Sections concerning the stationary regime and the departure from it occurring in 
the phase transition processes is centered on the behavior of the quantities tql(T), Cgl(^) an< ^ Vgl(T) defined in 
Eqs. (|Cip . These quantities, in turn, depend on the behavior of |a(T)|, which, as we have observed in Section Hill 
has the meaning of the mass term in the equations for the macroscopic wave function (the order parameter) er, or 
the wave function x normalized to |<t°|, namely x — l fJ l/l cr °l = {N/N ) 1 / 2 satisfying the general TDGL equation ((31))) 
(equivalent to Eq. ([33])). We might consider non- instantaneous phase transition processes where a time-dependent 
temperature might still be defined. In these processes, assume that the transition starts at the critical temperature 
Tc and the stable configuration is reached at the so-called "Ginzburg temperature" Tc, with Tq < Tc, after a certain 
interval of time during which the system is said to be in the critical regime [33| • In practice, the model here adopted 
aims at describing brain functioning as a repeated passage through a critical temperature state towards which the brain 
naturally relaxes after having reached a lower-temperature stable configuration, that corresponds to the response to 
some external stimulus or any other kind of perturbation. The crucial point is that this "critical dynamics" is always 
performed as a transition between different thermodynamic equilibrium states (at Tc and at Tq): irrespectively of 
the nonequilirium processes involved into this transition, this implies that the general thermodynamic description 
derived from the fluctuation theorem [34[ applies to this model of critical dynamics of the brain. The rate of entropy 
production that is inherent in the process described by the TDGL equation can be accordingly related to the amount 
of heat dissipated by the brain to reach new equilibrium configurations below the critical temperature Tq- 

We want to model time dependence of |a(T)| during such a critical regime evolution. 

Here we remark that the departure from the stationary regime (at Tc), representing the system ground state at 
a given time in the system evolution 2 , namely the start of the critical regime, is driven by fluctuations which may 
be spontaneous (typically occurring in a quantum vacuum) and which can trigger a phase transition, if the requisite 
transition energy is provided by some endogenous or external stimulus. These ground state fluctuations turn into 
temperature fluctuations since in the dissipative model the ground state is in fact a thermal state 0,13 [H. At the 
end of the critical regime (at Tc) the systems arrives at a "new" ground state configuration and the "phase transition" 
is thus completed. As a matter of fact, the brain system undergoes a continuous sequence of phase transitions in a 
path through the (infinitely) many ground states whose existence is allowed by QFT 0, Q [H[ . 

This is a delicate point, where the dissipative many-body model turns out to be very helpful in depicting what 
a traditional non-equilibrium thermodynamical approach cannot explain. It is well known that mammalian brains 
operate at a steady state constant temperature that is homeostatically regulated. Then it is crucial reconcile the 



2 See Refs. 0, 0, for the discussion of the existence of infinitely many vacua (ground states) representing different, physically 
incquivalent, microscopic configurations of the system. 
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temperature invariance of brains with the temperature fluctuations intrinsic to "far-from-equilibrium" activity seen 
in laboratory experiments as well as in everyday life. (As noted in Section I, brains consume free energy at rates 
ten-fold greater than any other organ caused by the heat exchanges at a local molecular level due to the ceaseless 
electrochemical and metabolic reactions). The resolution of the apparent contradiction between homeostasis and 
the far from the equilibrium brain activity is implicit in the dissipative many-body model, which suggests that a 
phase transition entails a brief, local fluctuation of temperature during the transition that punctuates the pseudo- 
equilibrium steady state. A peak is followed by reversal (a so-called biphasic transient), leading the system to a state 
that differs from the one before, as observed and predicted by the model. (For brevity we do not report more on such 
a hysteresis-like property in brain dynamics; see e.g. fA [35l - l38j V 

In fact, such a transient fluctuation in temperature has been demonstrated for nerve and modified muscle otherwise 
held in thermal steady state in experiments with the squid giant axon and siphon (39[. It shows that an action 
potential is accompanied by brief cooling of the nerve and synapse when the sodium ions expand into the intracellular 
compartment and the potassium ions expand outwardly. The cooling is soon overshadowed by the heat released by the 
burst of metabolism that pumps the ions back to their compartments and restores the energy reserve. The thermal 
reservoir of the water provides the buffering and averaging, such that the dissipative system can undergo substantial 
local fluctuations in temperature, which are masked in the global heat production. Thus we see another crucial 
role of water by continually holding temperature constant in the average, although locally (in space and in time) 
fluctuating. Our reference to the dependence of the system on temperature has to be understood as its dependence 
on such localized transients. 

It is interesting to note that in completely physically different contexts (such as inneuroscience, condensed matter 
physics, particle physics and cosmology) intensive theoretical and experimental research has shown that extended 
objects with non trivial topology (as for instance vortices) appear during the critical transition processes (see e.g. 
[IE Hill 53) and persist for varying time intervals thereafter. In the present Section we are indeed going to discuss 
the vortex solution to the TDGL equation. 

Due to the nonlinearity of the TDGL, numerical simulations have been associated often to the theoretical analysis 
(see [33| and Refs. quoted therein). Let us write Eq. as 



1 d X ^2 1 

\ ' 7 

GL 



= V Y , 

dq\j dt v x e 



n 2n c 2 f mv 
(1 - X ) - fGL 



X , (40) 



which admits as a solution the (quasi-)static vortex solution for = ^7§f ~ ( m sucn a case, indeed, it is 

recognized to be the vortex equation) [l!| [35|. The condensate produces the classical vortex envelope <r(r). Since 
(To is temperature dependent, the vortex envelope changes with T. The phase transition and the vortex formation 
appear thus as the effect of localized (non-homogeneous) boson condensation. 

Actually, the phase transition begins with an abrupt, adiabatic, localized decrease of the order parameter (a measure 
of the analytic power of the background activity) to near zero, denoted as the null spike, resulting locally in loss of 
order and coherence, hence onset of symmetry. Concomitant to this, the spatial variance of the analytic phase 
increases and a discontinuity in space and in time of the analytic phase appears. The start of the phase transition is 
thus adiabatic and it appears to be instantaneous at the time scale of the beta and gamma cycle durations. 

The extreme spatiotemporal localization of the null spikes indicates that they are associated to singularities. These 
constitute the core out of which, in the non-instantaneous phase transition process, the vortices start to form. The 
singularities coincide with the apex of the phase cones which start to develop and spread over the system. We can 
thus say that the null spikes mediate or precipitate the phase transition toward the formation of a new AM pattern. 
The non-homogeneous boson condensation at the origin of the formation of vortices and ordered localized patterns 
requires energy (ordered states are lower in energy and separated by an energy gap from the symmetric state). 
Consequently, apart the adiabatic instantaneous start, the phase transition itself appears to begin with cooling and is 
non-instantaneous. The cooling appears as a manifestation of the process of symmetry breakdown eventually leading 
to the emergence of new extended AM and PM patterns. It is thus different from the very localized cooling that may 
occur due to localized (e.g. on the axon) metabolic activity. 

It is interesting to observe that these features predicted by the model find a correspondence in the features observed 
at the level of neuronal assemblies. It appears that the neural transition requires three conditions as precursors. The 
first is macroscopic arousal with a high level of background 1/f activity (see our remarks in the following Section). 
Second is activation of a Hebbian nerve cell assembly, which amplifies, generalizes and abstracts a sensory signal into 
a burst of high intensity but microscopic firing of pulses. Third is the initiation of a narrow-band oscillation in that 
burst. These three changes require an increase in dissipation of free energy as heat and therefore in temperature. But 
the phase transition itself begins with cooling as predicted by the dissipative model. 

From the TDGL equation (see also Ref. [3| and [35]) one may also derive that, in the first approximation the vortex 
core size is of the order £q L cx |a| -1 . This means that the vortex core is increasing with temperature increase, as T 
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approaches Tc from below, the vortex envelope disappearing at Tc- Rising temperature above Tq leads to symmetry 
restoration (unstructured ground state). 

Vice-versa, going from above back to Tc the unstructured background activity (fully symmetric) with vanishing or 
very low analytic amplitude exhibits, at Tc, undefined analytic phase, namely the singularity (null spike) as it is at 
the center line of the vortex core. As T is lowered below Tc, the critical regime starts, vortices appear, whose core 
shrinks as temperature further decreases. 

As mentioned in Section II the singular NG boson (the phason) condensation is essential for obtaining a nontrivial 
topological charge. In Refs. (29|,|30( it is shown that the space-time-dependent order parameter can be gauged away by 
an appropriate gauge transformation when the boson condensation function is regular. When the boson condensation 
function presents singularities corresponding to regions occupied by the normal state (i.e. without condensation) the 
boson condensation is manifested in macroscopic structures. The vortex, e.g., is singular on the z-axis, at r = 0. In 
this case, the topological charge is characterized by the integer winding numbers n ^ 0. Associated to the vortex 
there is the quantized flux <j) = 2nn/q pi [l9l [35j. 

When the l.h.s of Eq. (|40]) is not neglected, one can show that in order to have a solution behaving as a oc 
exp — (tu/h^t, the condition 

k 2 > m 2 (t) , (41) 

has to hold during the critical regime between Tc and Tc for each fc-mode (fc = v^k 2 ), with k the wave number and 
to 2 (t) the time-dependent "effective mass" . In our discussion then we may follow the arguments presented in Ref. 
Q , where we have considered the formation of vortices during the critical regime and the appearance of null spikes 
(there we have considered the harmonic limit approximation consisting in neglecting the nonlinear term in the TDGL 
equation [11, E3, H3 ) ■ 

Suppose the critical regime starts and ends at the times t = and t = r, respectively. For a given k, Eq. (|4T|) holds 
up to a time Tfe such that m 2 (t) is larger than k 2 for t > r/.. The corresponding fc-mode can propagate in a span 
of time < t < Tfe, which determines the dimensions to which the domains can expand and the "effective causal 
horizon" [40L l4lj will be inside the system (possible formation of more than a domain) or outside (single domain 
formation) according to whether the time occurring to the fc-mode to reach the boundaries of the system is longer or 
shorter than the allowed propagation time, respectively. 

In order to determine the value of r^, one must assign the explicit form of m 2 (t). Its time-dependence may be fixed 
in such a way to allow vortex formation [TtJ • One possible modeling of m 2 (t) is then: 

m 2 (t) = to 2 e 2h ^ . (42) 

Let £ denote the correlation length corresponding to the fc-mode propagation and L oc m^ 1 . Then the correlation 
propagation time is implicitly obtained from Eqs. (|4"Tj) and (|42l) as: 

h(T k ) = In (—) oc In (^) . (43) 



m oJ \£. 

We see that L acts as an intrinsic infrared cut-off. Small values of fc are indeed excluded by the non-zero minimum 
value mo of m. Consequently, long wave-lengths are not allowed. This means that only domains of finite size can 
be obtained, which is what happens in the brain activity. The linear size of the domains over which the correlation 
length extends at the end of the critical regime is of the order of Afc oc m~ 1 (rfe). 

It remains to choose the explicit analytic expression for h(t) in Eq. 142[) . The choice we may adopt for h(t) is fl7j : 

where a,b,c are (positive) parameters chosen so to guarantee dimensionless h(t). Let A be an arbitrary constant and 
define the ratios c/a\ = tq, aX/b = t . We have H(tq) = h(ro) and the time derivative of h[t) (and thus of m 2 (t)) is 
zero at t = r = ±^tqTq, which thus plays the role of the equilibrium time scale. We have 

&(t) = ±J_-J__tw±Et ) (45) 
Xt Q 1 + 77 2 

for t 2 /t 2 ss 1 and T = 1/Xtq. In the linear approximation one finds that the number of vortices n^e/ possibly 
appearing during the critical regime is given by [171 141| : 

n def oc m 2 (r) w to„ \t/Xtq\ . (46) 
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Since the size of the vortex core is given by (m(£)) -1 , Eqs. (|42]) and (|45|) show that such a size evolves in time 
as e Tr *, t < t (t < Tfc for the fc-mode). Thus, as observed in Ref. |4|, we have both, converging (imploding) 
and diverging (exploding) regimes, which is in agreement with the phase cone behaviors observed in laboratory. 
The "normal" (disordered) state is confined to the vortex core, then in the case of enlargement of the vortex core 
(exploding regime) disorder (local correlations) prevail. Instead, the shrinking of such a region (imploding regime) 
may signal that ordering (long range correlation) is prevailing (the vortex is "squeezed out"). This is in agreement 
with laboratory observations which show that explosion or implosion is obtained if the local connections or the long 
axon connections predominate, respectively 42]. Observations also show that many phase cones exhibit little or 
no rotation but repetitive outward or inward pulsations with each cycle. When rotational gradients (vortices) are 
observed, the singularity is associated to the vortex core singularity. All four types of these observed spatiotemporal 
phase gradients (See Fig. 2) are thus predicted by the model. 

Finally, we note that conventional neurodynamics may explain the negative gradient (e.g. in terms of a pacemaker), 
but not the positive gradient. Moreover, in the conventional framework there is no explanation of why both gradients, 
the positive and the negative one, occur, one or the other at random. Furthermore, we offer as a prediction the 
expectation that the locations of the initiating null spike, the apex of the following phase cone, and the center of 
rotation or expansion-contraction of the vortices will be closely related though not necessarily coincident. At present 
the available experimental techniques indicate that this occurs, but they lack sufficient resolution and precision to 
provide experimental proof. 



V. POWER LAWS AND COHERENT STATES 



In this Section we consider one more of the predictions of the many-body dissipative model which is confirmed 
in laboratory observations. According to the model, the system ground state and the phase transition processes are 
characterized by the boson transformations, such as the one of the NG field described by Eq. (fT4"]) (see also Eq. ([T2"|l ). 
Boson transformations are the ones by which (Glauber) coherent states are constructed [IHESI- Coherence is therefore 
a distinctive feature of the basic dynamics of the system and one may expect that it manifests itself at a mesoscopic 
and macroscopic level since coherent states are known to be the quantum states most near to classical (macroscopic) 
states. Elsewhere 0, 0, 0] it has been stressed that one of the merits of the dissipative many-body model is that 
"classicality" emerges naturally out of the basic quantum dynamics, not as an artificial classical limit imposed by 
hand. The laboratory observations are in fact at a classical macroscopic level, as it should be (it has to be stressed 
that in our analysis neurons, glia cell and other biological structures are classical objects; the quantum variables are 
the dipole vibrational modes). However, some of the system features are macroscopic quantum features, in the sense 
that they cannot be explained without recurs to the underlying quantum dynamics. Among such features, there are 
the power-laws observed in the brain functional activity, on which we focus our attention in this Section. 

In the frame of the entire analytical function formalism, power-laws have been shown to be related to coherent 
states p4l [45l]. Let us very briefly summarize such a result. Power- laws denote scale free, self-similar phenomena of 
fractal nature characterized by a straight line graph in log-log coordinates, with the fractal dimension d corresponding 
to the slope of the line (the tangent y/x of the straight line angular coefficient). Self-similarity is a characterizing 
fractal property [HI, H?} expressed by 

(97)" — 1) f° r an y integer n, (47) 

where for convenience (see below) the number q is parameterized as q = l/f d and 7 is the parameter characterizing 
the specific considered fractal, whose dimension d is then given by 

(48) 

Log} 

This equation clearly describes a straight line in the log-log coordinates (x = Logf, y = Log"/). Here we denote logio 
by the symbol Log. We now remark that (97)™, up to the normalization factor l/vnl, is the nth element of the basis 
of the entire analytical functions defined in the space of the complex numbers z = 37. In Eq. (1471) the number qj is 
considered to belong to the restriction to the real axis of the z-complex plane. The connection with coherent states is 
then established by observing that they are in fact constructed in the Fock-Bargmann representation (FBR) by use 
of the entire analytical functions [19j, |43( . More precisely, the coherent states for such a z = q 7 variable are defined to 
be "q-deformed" coherent states of coherence strength 7, or, also, squeezed coherent states, q being the deformation 
or squeezing parameter. We will not insist on such technical points here and the reader may consult Refs. [HI, SH| for 
details. The point is that, as a result, self-similar processes of fractal dimension d can be associated to g-deformed 
quantum coherent states, namely they appear as macroscopic systems generated by local deformations of coherent 
quantum dynamics, they are macroscopic quantum processes induced by quantum boson transformations. 
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FIG. 3: A. Multichannel recording from a high-density array fixed on the auditory cortex of a rabbit at rest. B. The power 
spectral densities in the working ECoG were computed for all 64 signals and averaged. The power-law trend lines (1//° and 
I// 2 6 ) were drawn by hand to emphasize the multiple peaks of power in the theta and beta-gamma ranges above the line, 
which were missing in the resting ECoG. From Figs. A1.01 and A1.02 in [13] 



Let us now turn to the laboratory observations (see e.g. |9^ll2l. |48||) and to the graphs of the power laws depicted 
in the lo-log coordinates in Figure [3] and |U On the basis of what said above, it is evident that we face macroscopic 
quantum processes characterized by a coherent dynamics. Coherence is in fact the property to which the analysis in 
terms of neuronal correlates also leads. This can be seen as follows. In laboratory observations [9T-fl2l l48j multichannel 
recording yields a collection of state variables. Weakly coherent signals are obtained in ECoG with electrodes which 
are sufficiently far apart. The cortex may be considered to be modular, and each module gives a signal. Coherent 
signals are instead obtained when the electrodes are closely spaced. In such a case we may talk of variables defining 
a state space. They are all from the same module, and their variations in time define a trajectory through the 
modular state space. Many data have been collected [9Ttl2| on ECoG spatial imaging coming from small, high-density 
electrode arrays fixed on the surfaces of the olfactory, visual, auditory, or somatic modules (typically 8x8 electrodes, 
spaced 0.79 mm apart, giving a 5.6 x 5.6 mm square aperture onto the ECoG). Figure [3] A shows the unfiltered 
background ECoG, obtained in such a way, which is featureless in appearance. Spectral analysis of the ECoG shows 
broad distribution of the frequency components in Figure [UB, where the temporal power spectral density (PSDt) 
in log- log coordinates is power-law, 1// Q , in segments. Below an inflection in the theta-alpha range the PSDt is 
flat, a = 0. Above, the logio power decreases linearly with increasing logio frequency in the beta-gamma range 
(12.5 — 80 Hz) with the exponent a between 2 and 3. One can show that in slow wave sleep the exponent averages 
near 3 [4{|; in seizures it averages near 4. Above 75 Hz the slope either increases (a = 4, 5{|), or it flattens to a = 0. 
On the basis of what said above, such values of the slope a provide corresponding values of the fractal dimension 
d = a, with deformation parameter q = 1/ f a and coherent strength 7 corresponding to the power spectral density. 
In Figure (3JB, each peak of power above the 1/f trend line in the beta-gamma range (arrows in Figure OB) reflects 
a brief epoch of narrow band oscillation. These multiple peaks indicate a departure from the scale free regime (the 
straight line) and therefore the presence of structures emergent from the background activity. Figure HjA, B shows 
the post stimulus time histograms (PSTH) of the microscopic pulses from a representative neuron in an excitatory 
population due to low and high intensity electric shocks. The averaging gives the impulse response of the macroscopic 
population. It presents a rapid rise and an exponential decay to the background firing rate. When the shock strength 
is reduced to threshold, the decay rate approaches zero. In dynamical terms the asymptotic convergence of activity to 
the background after perturbation is evidence that a point attractor governs the cortical background activity at unity 
gain (pp. 285-305 in |48j). An attractor by definition represents a stable (ground) state. The power law in Figure 
IHC, D signals the scale free ("spatial similarity") property of such ground states. The spatial similarities reveal the 
long spatial distances across which synaptic interactions can sustain coherence, namely the high-density coordination 
of neuron firing by synaptic interaction, in agreement with the dissipative model prediction. In the frame of our 
discussion of theoretical and observational results, the coherent mutual excitation among pyramidal cells thus appears 
as an especially attractive conclusion, also considering that they comprise about 80% of cortical neurons, and they 
provide 90% of the synapses of cortical origin to each other. 
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FIG. 4: A, B. The symbols A show the post stimulus time histogram (PSTH) of a representative neuron in response to a 
weak electric shock, fitted with the solution to a 4th order linear differential equation [4g|. The prolonged discharge without 
inhibitory overshoot is reverberation that is due to mutual excitation. C, D. The decay rate determines the inflection frequency 
of the power-law PSDt, and the rise rate determines the exponent, a. The predicted range of the trend lines is 2 < a in 1/ f a . 
We propose that the segments give the canonical form of resting cortical PSDt- From Figs. 5, 7 in 



VI. CONCLUSIONS 



The brain is an open system submitted to continuous exchange of energy and information with the environment. 
Sensory inputs drive the brain out of its background activity and large amplitude modulated patterns appear [§-[l2j 
which are the expression of long range neuronal correlation. Moreover the brain manifests a form of conditional 
stability, criticality, denoting a state of readiness to transit from one state to another, like the one of a subject who is 
expecting one of several CS, by which to choose one of several courses of action. In the expectant state of searching, a 
sensory cortex can be viewed as having a set of attractors, one for each expected CS, and each with its basin around 
an attractor in an attractor landscape [511 ] . The CS determines the choice by directing the state trajectory into a 
particular basin. Convergence to the attractor directs the cortex to form the AM pattern that is selected by the CS. 
The cortical assemblies of neurons can complete widespread phase transitions in very few milliseconds, regardless of 
their correlation lengths and carrier frequencies. This property can help explain the rapidity of the flow of mental 
images. The distances across which coherent oscillations are maintained (the correlation lengths) are far larger than 
the modal diameters of the dendrites of the participating neurons and can include multiple cortical areas, even the 
entire scalp, which can help to explain the multisensory integration that is required for Gestalt formation. Moreover, 
spectral analysis of the ECoG shows the occurrence of l// a power-law form of the temporal spectral density in log- log 
coordinates which indicates that the activity is scale-free. 

The above properties are the macroscopic manifestation of the mesoscopic and microscopic dynamics. The formation 
of AM and PM assemblies of neurons have been associated, in previous works [3, @] , with the spontaneous breakdown of 
symmetry occurring at the level of the many-body system components and their coexistence and temporal succession 
to the dissipative character of the dynamics. The large distances across which coherent oscillations are observed to 
occur in the brain is accounted by the long range of the correlation, extending even to the whole system volume, 
predicted as a consequence of the spontaneous breakdown of symmetry in QFT. In this paper, continuing our study 
of the many-body dynamics underlying these behaviors, we have focused our attention on the non-equilibrium non- 
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instantaneous phase transitions occurring at the level of the many-body dynamics. In the stationary regime (Section 
HI)) , as customary, one considers the extremizing condition of the system free energy J^- = whose solutions describe 
the ground state of the system. In the brain activity, however, the brain is continuously moved out from its ground 
state activity entering a non-stationary dynamical regime. We have thus considered (Section IIIip a non-vanishing 
value for the quantity J^- , which now expresses the rate at which the system approaches to the stationary regime. 
In such a way we have obtained the time-dependent Ginzburg-Landau equation. During such a non-equilibrium 
phase transition, named the critical Ginzburg-Landau regime, the system dynamics turns out to be characterized by 
topologically non-trivial structures (vortices) , which in the formalism of the dissipative model of brain are described by 
non-homogeneous boson condensation processes. The size and the life-time of these topological structures have been 
discussed (Section llVp and their behavior commented upon in relation to the temperature changes in the domain of the 
critical regime, from the critical temperature Tq at which the transition starts to the so-called Ginzburg temperature 
Tq, with Tq < Tq- The number of vortices appearing in the phase transition process has been estimated, in agreement 
with a previous analysis In such a critical phase transition the system moves towards a new stationary ground 
state, which thus describes an attractor for the system non-linear evolution. Since our frame is the one of QFT, which 
posses infinitely many unitarily inequivalent ground states, the critical system transition to the stationary re gim e is 
not at all described by a trivially determined trajectory in the state space of the system. Elsewhere [3, 0, [36- 38] we 
have shown that such a trajectory is indeed a deterministic chaotic trajectory in an attractor landscape, the initial 
conditions being determined by the CS specificity. 

Finally, the observation of the l// a power-law of the temporal spectral density in the lo-log coordinate plot, which 
exhibits a scale free activity in the brain ground state activity, has been described, in agreement to previous studies 
of the relation between fractals and coherent states Ell, as the macroscopic manifestation of the underlying 
coherent many-body dynamics, namely in terms of the scale-free property of coherent states. In such a derivation, 
the properties of the entire analytical functions have been used. 



Appendix A 

The ECoG and EEG present nearly random but bounded amplitude and spectral distributions, often power-law 
(1//). The oscillations in the beta-gamma range occur in brief epochs with spatial patterns of amplitude and phase 
that are statistically constant. In those stationary epochs, the two major operations of normal cortical dynamics of 
dendritic integration of waves and axonal transmission of pulses are executed in near-linear domains (48j . Superposition 
justifies our use of the tools of linear analysis. The interactions among excitatory neurons that produce the random 
1// oscillations in the ECoG and EEG background activity at rest has been modeled with a single positive feedback 
loop [49l | . and the interactions among excitator y a nd inhibitory neurons that produce the beta-gamma oscillations 
(12 — 80 Hz) with a negative feedback loop [48l.]52T|. 

Multiple AM/PM patterns coexist in different frequency bands, overlapping in space and time. Whether they 
interact and if so how strongly has not been determined; they appear to have something like the imperturbability of 
solitons. The stability of the patterns embedded in and arising from the background turbulence leads us to view cortical 
dynamics as bistable. In the behavioral state of expectancy the cortex is in a receiving mode that is characterized 
by disorganized, random pulse firing and chaotic wave activity. We adopt the AM pattern feature vector as an order 
parameter. For such a state it is zero, and say that the cortical populations are in a gas-like phase. In the behavioral 
state of perceiving the cortex is in a transmitting mode of strong interaction, which imposes order on the firing of 
pulses and synchronization on the carrier wave. The order parameter is non- vanishing positive. The change resembles 
the phase transition from a gas-like phase to a liquid-like phase of condensation. Usually on close examination of the 
analytic phase we can find a discontinuity in the time series at which the new AM/PM patterns begin Q suggesting 
a Type I phase transition. After ^3 — 5 cycles of the carrier the patterns more slowly dissipate, suggesting a Type 
II phase transition. 

The transitions only occur when a CS is presented and the subject is in an expectant state. We believe that the 
phase transition requires a non- vanishing activation energy, strong enough in order to minimize the risk of false alarms 
and to conserve energy, and that this energy is provided by a Hebbian assembly, which is an intracortical mesoscopic 
network that has been formed by prior reinforcement learning each discriminable CS. Its functions are to amplify, 
abstract and generalize each weak sensory input by reverberation (mutual excitation) to activate the entire asse mbly 
by any part. The Hebbian assembly also greatly facilitates narrow band oscillation in the beta-gamma range [48j . 
The bandwidth is not zero; the narrow distribution about the center carrier frequency causes beats, which appear as 
reductions in amplitude that recur at intervals proportional to the bandwidth regardless of the center frequency [53j . 
During a beat the power in the designated frequency band can briefly go to zero [Hi| in one or more null spikes j4j. 
We propose that the null spike manifests a space-time singularity at which the phase transition begins. Its location 
may weakly correspond to those of the apex of the phase cone and the center of vortex rotation. 
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Evidence for the vector fields is derived from recording and interpreting the highly textured spatial patterns of 
oscillations provided by the images of scalar fields of potential from the EEG and ECoG. The neural mechanisms 
of vector field formation are modeled with large networks of differential equations that we solve by piece- wise linear 
approximation (Chapter Six in [48[ and by implementation in VLSI 55]. The solutions show the existence of a 
point attractor that homeostatically regulates the background activity that maintains cortex in a stable state of 
criticality. This is the receiving phase [49( , in which the oscillatory impulse response has en exponential envelope with 
an exponent, a < 1. At criticality a = 0, endowing cortex with exquisite sensitivity to selected input, form scale-free, 
power-law distributions of its parameters and state variables, and establish remarkably long correlation distances that 
can include the entire cortex of each hemisphere. In the transmitting state a > 0, and the cortex transits to a limit 
cycle attractor in a form of subcritical Hopf bifurcation. It is the amplified output of the Hebbian assembly that 
selects the carrier frequency, and it is the null spike that drives the cortex inexorably to the limit cycle attractor. 



Appendix B 

Let ip(r, t) — \J a(r, t) exp(iS/h) denote the wave function for the system quantum components, i.e. the microscopic 
wave function. It satisfies the Schrodinger equation from which the continuity equation 

| + V'H = (Bl) 

can be derived. As well known, by defining 

1 h 2 vV/ 2 

v=-VS and Q = _IL-5L^_, (B2) 
one obtains the hydrodynamic equation of motion (the Madelung picture) 

^ + (v- V)v = --V(F + Q), (B3) 
at m 

where V is the potential appearing in the Schrodinger equation and Q is the so-called quantum potential. In the 
presence of a vector potential A the (hydrodynamic) velocity v is given by 

v = — CVS-qA). (B4) 
m 

We see that the chemical potential fjb\ in Eq. (|2ip is equivalent to the quantum potential Q in Eq. (|B2|) . 
We consider now the derivation of Eq. (|22l) . We apply the operator V to both members of Eq. ([20]) and use the 
identities: 

2 

TftV 

V—— = m(v • V)v + mv x (V x v) , (B5) 



dv dv , ov v „ 

_ s _ + (v.V)v= w + V T -vx(Vxv). (B6) 



where v = v • v. The rotor of Eq. (|B4I) gives 



V x v = — — V x A = — — B. (B7) 



Combining these relations we get Eq. 



By adopting the notation of Ref. [3l| , we put 



Appendix C 



TGL(T) = R7T 5gl(T) = 2m\a(T)\' mh{T) = Rl^ = iV°(T) ' (C1) 
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and using Eq. ([26|) and 
Eq. (gT]) takes the form: 

Use of the relations 



R ee £| L ( V + JrA^-ryGLlcf + l, 



. , Oct - .. R 

in—- = Ha + in a. 

at T G L 



pl/2 7VV2 

then leads us to Eqs. (JSgJ and O- 

The (stationary) case tgl — > (i.e. 7 — » 0), can be studied by considering Eq. (|27l) . Multiplying both members 
times 7, we see that the limit 7 — ?> leads to the stationary GL equation (|2~5l) . On the other hand, by multiplying 
both members of Eq. (|C3|) times tgl , the limit tgl — > leads to 

2 

a = (C5) 



V-v = 



(C2) 



(C3) 



(C4) 



Ra=-[G- i&l 



V J 



ft p 



The real part of Eq. (|C5|) gives G cr — 0, i.e. 



'-&(x)' 



X - X - 0, 



where \ = l cr l/l (j0 | = (N/N ) 1 / 2 . The vanishing of the imaginary part of Eq. (|C5|) implies 

V-J = 0, i.e. v-ViV = 0, 
which means that in the stationary case non-homogeneities in N only arise in the plane perpendicular to v. 



(C6) 
(C7) 



Appendix D 

We can define fi = fio + Ml + M2 and introduce the total chemical potential (cf. Eqs. ([29]) and (|30|) ) 

mi) 2 / ft 2 V 2 ^ 1 / 2 £ 2 m V • J\ tow 2 , , 

Me =M+^- = Mo-^ ^T7^-— + ^7~' D1 

2 \ 2m iv 1 /^ tgl q N J 2 

The gauge invariant potential $ is defined as 

86 ( mv 2 \ 
$ = < w , + ft_ = -U+_J=- A , B . (D2) 

and we see that /x B is always opposite in phase to $. We may also introduce £ = fi + q(p and the electrochemical 
potential £s given by 

tow 2 t <9# , T ,„. 

Cs = C + — = Ms+9<^= -ft^. (D3) 

We see that it depends on the rate of change of the NG mode (the phason mode 9). Finally, the gauge invariant 
vector and scalar potentials are given by 

tov $ 
A g = — and cpg = — j (D4) 

respectively. The corresponding electric field is 

E = -^-V^ g (D5) 

The Coulomb gauge condition V ■ A g = then implies that 

V • v = (D6) 



and Eq. (IB7|) is also reobtained. Due to Eq. (|D6|) . the velocity normal to any boundary thus vanishes, n ■ v = 0, and 
at boundaries any current such that n-v/0 converts into a current with n ■ v = and the condensate density N 
has nonuniform distribution. 
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